Electron correlation in FeSe superconductor studied by bulk-sensitive photoemission 

spectroscopy 
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We have investigated the electronic structures of recently discovered superconductor FeSe by 
soft-x-ray and hard-x-ray photoemission spectroscopy with high bulk sensitivity. The large Fe 3d 
spectral weight is located in the vicinity of the Fermi level (Ef), which is demonstrated to be a 
coherent quasi-particle peak. Compared with the results of the band structure calculation with 
local-density approximation, Fe 3d band narrowing and the energy shift of the band toward Ef are 
found, suggesting an importance of the electron correlation effect in FeSe. The self energy correction 
provides the larger mass enhancement value (Z~ x ~3.6) than in Fe-As superconductors and enables 
us to separate a incoherent part from the spectrum. These features are quite consistent with the 
results of recent dynamical mean-field calculations, in which the incoherent part is attributed to the 
lower Hubbard band. 



PACS numbers: 79.60.-i, 74.25.Jb, 74.70.-b, 71. 20. Be 

I. INTRODUCTION 

Fe-based high-T c superconductors have attracted enor- 
mous attention for their possibly new-type superconduct- 
ing mechanism and the potential of breaking the deadlock 
in the high-T c superconductor research field. A fluorine 
doped LaFeAsO has been discovered to be a supercon- 
ductor below T c — 26 K, which contains the two dimen- 
sional Fe plane in the Fe2As2 layer pi So far, several tens 
of superconductors in limited types of mother materials 
such as LaFeAsO, BaFe2As2, and LiFeAs, have been syn- 
thesized^ In addition, a simple Fe-Se binary compound 
(FeSeo.82) has been discovered to show superconductiv- 
ity^ The appearance of superconductivity in the FeSe 
system indicates Fe2A2 (X=P, As, and Se) layer is es- 
sential for the superconductivity in these Fe-based su- 
perconductors. The density functional study has pointed 
out that FeSe is not a conventional electron-phonon su- 
perconductor, being similar to LaFeAsOi-zF-r system* 4 . 
Other common features to the Fe-based superconductors 
have also been revealed, which both the antiferromag- 
netic spin fluctuation and the anion height are closely 
related to the appearance of the superconductivity^ 

There are, however, many differences between FeSe 
and other Fe-based superconductors: in FeSe (i) there 
is no separating layer and the Fe2Se2 layer is electri- 
cally neutral, (ii) the superconductivity is very sensitive 
to the deviation from the stoichiometric composition* 7 ^ 



(iii) there is no magnetically ordered state in p-T phase 
diagram. 9 Furthermore, it is pointed out by theoreti- 
cal studies that the electron correlation effect in FeSe 
is stronger than other Fe-based superconductors ! 10 ' 11 In 
FeSe, only a few experimental results, for instance, Som- 
merfeld coefficient 7 (= 5.4-9.1 mJ/mol K 2 ) have been 
so far reported 3 -' 7 - although the electron correlation ef- 
fect in other Fe-based superconductors was investigated 
by spectroscopic experiments in detailJ^— . The reason 
why no angle-resolved photoemission spectroscopy exper- 
iment has been reported to quantitatively evaluate the 
electron correlation effect in FeSe is that a high quality 
single crystal is difficult to be grown. Recently, Fe(Se,Te) 
system has been intensively investigated due to the suc- 
cess of the high-quality single crystal growth. Even so, 
the Fe(Se, Te) system is essentially different from the end 
member FeSe in the sense that the magnetically ordered 
state and remarkably large 7 value (= 39 mJ/mol K 2 ) 
have been foundj 16 ' 17 

We have examined two different synthesis processes 
for the FeSe superconductor and performed the soft-x- 
ray and hard- x-ray photoemission spectroscopy (SXPES 
and HAXPES) in order to quantitatively evaluate the 
electron correlation effect in the bulk. The SXPES and 
HAXPES have been widely recognized as the powerful 
techniques which can reveal bulk electronic structures 
due to the long inelastic mean free path of photoelec- 
trons excited by high-energy x rayJ^— It is found in the 
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angle-integrated PES spectrum that a large Fe 3c? spec- 
tral weight is located in the vicinity of the Fermi level 
(Ep) and it decreases steeply toward Ep, being a simi- 
lar feature to those in the other Fe-based superconduc- 
tors non-doped LaFePO and LaFePnO0.94F0.06 (Pn=P, 
As) i 12 ' 23 Considering the self energy correction to the re- 
sults of band structure calculations, the experimentally 
observed band narrowing and the energy shift of the band 
toward Ep are fully explained. The correction also pro- 
vides the renormalization factor Z of ~0.28 and enables 
us to separate the incoherent part of the quasi-particle 
spectrum. 



II. EXPERIMENTAL 

For SXPES and HAXPES measurements, single- 
crystalline tetragonal [Tet(Single)] FeSe and polycrys- 
talline tetragonal [Tet(Poly)] FeSe were employed, to- 
gether with polycrystalline hexagonal [Hex(Poly)] FeSe 
(likely, FerSes) as a reference material* 7 - The Tct(Singlc) 
FeSe was grown by the chemical vapor transport method 
using Fe and Se powders^ It was found to contain 
the hexagonal phase around the crystal edge by x-ray 
powder diffraction (XRD) measurements. The coexis- 
tence of two stable phases has also been reported by 
other groups.— i2£r— The spot size of the SX beam in 
PES measurements is small (~ 10/im x 100/im) enough 
to exclude signals from the hexagonal phase near the 
crystal edge because of the much larger sample size 
(~ 500/xm x 500^m). The Tet(Poly) FeSe was synthe- 
sized using a high-frequency induction furnace. The fur- 
nace is very useful for quick cooling of the samples be- 
cause this furnace can heat only the samples. This makes 
the pass time through other low temperature phase rel- 
atively short. The Hex(Poly) FeSe was prepared by a 
conventional solid state reaction. Crystal structures of 
both the Tet(Poly) and Hex(Poly) FeSe were also eval- 
uated by XRD. The Rietveld analysis against the XRD 
profiles reveals that the Tet(Poly) FeSe sample contains 
the hexagonal-phase FeSe of 10% as an impurity and the 
Hex(Poly) FeSe sample has a single phase with in the lim- 
its of the resolution. The Tet(Single) FeSe and Tet(Poly) 
FeSe have the transition temperature T c zoro ~6 K esti- 
mated by p— T and T c ~8 K estimated by \ — T measure- 
ments under ambient pressure, respectively, which are 
similar to reported values! 3 ' 7 ' 25 This implies the present 
samples have the selenium defect of few percent. 

SXPES was carried out at both the Figure-8 undula- 
tor SX beamline BL27SU in SPring-8 using the SPECS 
PHOIBOS 150 hemispherical electron energy analyzer 
and the twin-helical undulator SX beamline BL25SU in 
SPring-8 using GAMMADATA-SCIENTA SES-200 spec- 
trometer^^ The highest total energy resolution AE 
[the full width at half maximum (FWHM)] was set to 
75 meV at hi/=600 eV. HAXPES was performed at the 
beamline BL19LXU in SPring-8 with MB Scientific Al- 
HE spectrometer. The linearly polarized light was de- 
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FIG. 1: (Color online) Se 3d core-level PES spectra of FeSe. 
Circles and thick solid curve indicate the experimental and 
fitted spectra, respectively. The Se 3d components of tetrag- 
onal and hexagonal phases are shown by thin solid and broken 
curves. The other components, for instance, Fe 3p and back- 
ground are not displayed for simplicity. 

livered from an in- vacuum 27 m long undulator. 30 The 
AE for the valence-band PES spectrum near Ep was set 
to ^380 meV at \iv ~8 keV. Clean surfaces of the sam- 
ples were obtained by fracturing samples in situ in UHV 
(^4xl0~ 8 Pa) at the measuring temperature (T ~20 K). 

III. RESULTS AND DISCUSSION 

FeSe has two stable crystal phases, that is, the tetrag- 
onal and the hexagonal phases and often contains both 
in one sample as has been mentioned above. In order 
to investigate the crystal phase mixing in samples, Se 3d 
(including Fe 3p) core-level PES spectra were measured 
for Tet(Single), Tet(Poly), and Hex(Poly) FeSe samples. 
The Se 3c? core-level PES spectra in both crystal phases 
excited by SX have a sharp doublet peak structure origi- 
nating from Se 3c? 5 / 2 and 3g? 3 / 2 components as shown in 
Fig. [TJ There is no surface-derived peak and only weak 
plasmon satellites are seen on the higher binding energy 
(Eb) side of the Se 3d peaks, indicating that the sur- 
face contributions to the spectra can be ruled out in the 
SXPES for the shallow core levels (including the valence 
band) for FeSe. In addition, the spectra in the tetragonal 
FeSe have the shoulder structure on the lower Eb side, 
which originates from the Fe 3p states. In the HAXPES 
the intensity of the shoulder is enhanced since the pho- 
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FIG. 2: (Color online) Valence-band PES spectra of FeSe. (a) 
Overall valence-band SXPES and HAXPES spectra. These 
spectra are normalized by the area under the curves after sub- 
tracting the Shirley-type background. 33 (The dotted line indi- 
cates the extrapolated one assuming the Lorentzian line shape 
in order to normalize the HAXPES spectrum.) (b) High- 
resolution PES spectra near Ef- The spectra of polycrys- 
talline samples are normalized so that the spectral intensity 
agrees with the intensity in each normalized high-resolution 
overall valence-band PES spectrum (shown partly by dots) . 



toionization cross section a of Fe 3p states relative to that 
of Se 3c? states increases by a factor of 8 compared with 
the SXPES.— We have obtained the tetragonal and/or 
hexagonal components in each sample by a deconvolu- 
tion procedure*^ The Se 3d peaks in the tetragonal and 
hexagonal components are located at the certain ener- 
gies in these samples, respectively, as shown by solid and 
dotted bars in Fig. [T] The peak of the tetragonal compo- 
nent is located at the Eb which is about 300 meV higher 
than that of the hexagonal one. This shift is caused by 
the structural difference, in other words, the difference 
of the covalent bond strength between these compounds, 
being consistent with what Se 4p states have a different 
structure in the valence band of these compounds (dis- 
cussed later). We note that the Tet (Poly) FeSe has the 
hexagonal component of 20 % in the Se 3d core-level PES 
spectrum. The influence of the hexagonal FeSe inclusion 
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FIG. 3: (Color online) Fe 2p core-level SXPES and HAXPES 
spectra of FeSe. 



on other spectra of Tet (Poly) FeSe will be discussed later. 

Figure[2](a) shows valence-band SXPES and HAXPES 
spectra of the tetragonal and hexagonal FeSe. In both the 
tetragonal and hexagonal FeSe the SXPES spectra have 
a sharp peak in the vicinity of Ep. In addition, there are 
some broad peaks and hump structures in the higher Eb 
region. The spectra of the tetragonal FeSe have the fea- 
tures very similar to the recently reported result.— These 
structures consist of Fe 3d and Se ip states as has been 
revealed later by comparison with theoretical spectra. In 
the HAXPES spectrum, it is found that the peak near Ep 
is suppressed and the structures in the Eb=3-7 eV be- 
come dominant. Such a variation of the spectral shape in 
the SXPES and HAXPES does not mainly originate from 
the increase of the bulk sensitivity but from the ratio of 
a between Fe 3d and Se 4p states, that is, Fe(3d)/Se(4p) 
is 2.8 at hv ~690 eV and 0.033 at ~8 keV. This indi- 
cates that the HAXPES spectrum of the valence band is 
nearly equivalent to the Se 4p (including weak and broad 
s) states. 

Now we focus on the similarity and difference of the 
electronic structures between the tetragonal and hexag- 
onal FeSe. They have rather different Se Ap electronic 
structures between Eb= 3 eV and 5 eV as recognized in 
Fig. [2] (a). Furthermore, high-resolution PES spectra re- 
veal that Fe 3d states near Ep also have different features 
as shown in Fig. [2] (b). We note that each spectrum of 
Tet (Poly) and Hex(Poly) FeSe is normalized by the inte- 
grated intensity of its high-resolution whole valence-band 
PES spectrum after subtracting the background. The 
Tet(Poly) FeSe has a prominent peak at Eb — 180 meV. 
The photoemission intensity decreases steeply toward Ep 
and is very weak at Ep. Meanwhile, the Hex(Poly) FeSe 
has a less prominent peak, which is closer to Ep. In 
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FIG. 4: (Color online) Valence-band PES spectra of FeSe 
measured at two different photon energies. The spectrum 
indicated by dots is difference between the off-resonance 
(measured at h^=690 eV) and anti-resonance (measured at 
hi/=702 eV) spectra for (a) tetragonal FeSe and (b) hexagonal 
FeSe. In the middle of Figure (a), calculated Fe 3d spectrum 
including the self energy correction and the broadened Fe 3d 
PDOS are shown by solid lines. The original DOSs calculated 
by Singh el al. are also shown.— 



addition, the obvious Fermi cut off is observed in the 
Hex(Poly) FeSe. Here we would stress that the spectra 
near Ep in the Tet(Single) FeSe is in good agreement 
with in the Tet(Poly) FeSe, indicating that the influence 
of the hexagonal FeSe inclusion in Tet(Poly) FeSe is neg- 
ligible near Ep. 

Fe 2p core-level PES spectra of the tetragonal and 
hexagonal FeSe are shown in Fig. [3] There is no charge- 
transfer satellite, as has been pointed out in LaFeAsO. 23 
In addition, one can see the sharp-peak-and-shoulder 
structure in the 2p 3 / 2 component of the SXPES spec- 
tra of both the Tet (Single) and Tet(Poly) FeSe as indi- 
cated by solid and dotted bars. Considering that the 
sharp peak becomes more prominent in the HAXPES 
spectrum, the peak can be assigned to the bulk \2p 5 3d 6 ) 
component. The sharp peak structure suggests the exis- 
tence of the highly coherent valence electrons in the bulk. 
Then, the shoulder structure is attributed to the photoe- 



FIG. 5: (Color online) High-resolution PES spectra near Ef- 
(a) The broken lines indicate the calculated Fe 3d spectrum 
(blue) and Fe 3d PDOS (black) without broadening and finite 
temperature effects. Other explanations are the same as in 
FigSl(a). (b) The PES spectrum measured at h^=40 eV. The 
calculated Fe3d spectrum without broadening is also shown. 

mission from the surface (and grain boundaries in the 
polycrystalline sample). These indicate that the photo- 
electrons with the low Ek (of less than a few hundred 
eV) lead to the large surface contribution to the spec- 
trum even in the SXPES experiment^ The possibility 
cannot be ruled out that the slightly stronger shoulder 
intensity in the SXPES spectrum of Tet (Poly) FeSe than 
that of Tet (Single) FeSe originates from the influence of 
the hexagonal FeSe inclusion. Meanwhile, the Hex(Poly) 
FeSe has a broad peak originating from multiplet struc- 
tures. The overall spectral shape is very similar to the 
reported spectrum in the hexagonal FeySeg^ 

Let us discuss Fe 3d states in the valence band PES 
spectra of the tetragonal and hexagonal FeSe. In Fig. @] 
(a), the valence-band PES spectra of the Tet(Poly) 
FeSe measured at two different photon energies are 
shown. The solid line indicates the spectrum obtained at 
hi/=690 eV, which is labeled as "off-resonance" in com- 
parison to another spectrum. The dotted line shows the 
spectrum which was measured at the energy just below 
the threshold of Fe 2p-3d absorption edge, labeled as "Fe 
2p anti-resonance" . In this spectrum Fe 3d states are 
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strongly suppressed since the tuned photon energy cor- 
responds to the energy providing the mostly minimum 
transition probability in the Fano lineshape.— One can 
see that the spectral weight between Ep and Eb—3 eV 
is remarkably reduced. This indicates Fe 3d dominant 
states are located in this Eb range. However, some spec- 
tral weights still remain, suggesting that there are Se Ap 
states hybridized with the Fe 3d states as seen in the 
HAXPES spectrum of Fig. [2] (a). Meanwhile, no signif- 
icant reduction above Eg=3 eV is seen except for the 
slight decrease around E B =A eV, since these structures 
mainly originate from Se Ap states. We note that one 
often employs the (on-) resonant PES to investigate the 
contribution of the specific electronic states by using the 
photon energy tuned to the core-level absorption maxi- 
mum. In the present case, however, the Auger decay pro- 
cess becomes dominant and the valence-band structures 
are smeared out due to the large background. The dif- 
ference spectrum between the off-resonance (hz/=690 eV) 
and anti-resonance spectra is also shown as dots in Fig.0] 
(a), which represents the Fe 3d partial density of states 
(PDOS) if the electron correlation effect is negligible. It 
has a sharp peak in the vicinity of Ep and broad hump 
structures at around Eb=1-5 eV and 4 eV. The differ- 
ence spectrum has similar features to the experimentally 
obtained Fe 3d states of LaFePnO0.94F0.06 {Pn=P, As). 23 
Meanwhile, the Hex(Poly) FeSe has much different Fe 3d 
states from the Tet(Poly) FeSe as shown in Fig. [4] (b). 
The difference spectrum has a large spectral weight up 
to Eb=6 eV and broad hump structure around 3 eV. 
These features are consistent with those in FeySes spec- 
trum measured at hi/=100 eVM- The SX spectrum has a 
stronger peak near Ep than the lower-hz/ spectrum pre- 
sumably due to the high bulk sensitivity. 19 

The results of the band structure calculations with 
the local-density approximation (LDA) for the tetrag- 
onal FeSe done by Singh et al. and the Fe 3d PDOS 
broadened by Gaussian and Lorentzian functions repre- 
senting the experimental energy resolution and the life- 
time broadening effect are also shown in the Fig. @] (a). 
By comparison with them, narrowing the Fe 3d band 
width and the energy shift of the band toward Ep are 
found in the difference spectrum of the Tet(Poly) FeSe, 
suggesting that the electron correlation effect cannot be 
negligible to discuss the electronic structures in this sys- 
tem. In order to take into account the electron correla- 
tion effect quantitatively, the self energy £(k, u>) defined 
as S(k,uj) = G ~ 1 (fc,w) — G^ 1 (k,ui) is often employed, 
where Go(k,u>) and G(k,uj) are the one-electron Green's 
functions without and with the electron-electron inter- 
action. We have calculated the spectrum with the self 
energy correction in accordance with the procedure re- 
ported by Shimada et al££- The fc-independent self en- 
ergy is assumed to be given by S(uj) — guj/(uj + iT) 2 , 
where g and P are employed as fitting parameters and 
huj < fi(= 0) for the occupied state. The fe-integrated 
one-particle spectral function for Fe 3d states including 
the electron-electron interaction, hereafter called "Fe 3d 



spectrum" Ad(uj), is obtained as follows, 
A d (uj) = ^2A d {k,uj) 

V / 

IT 



(1) 



de Q k D d (e° k )ImG(k^) 



de° k 



D d (el) 
lmS(k,uj) 



{Two - e° - Rei7(fc,w)} 2 + {Imi7(A;,a;)} 2 

where D d (e k ) denotes Fe 3d PDOS. The obtained spec- 
trum after optimizing the parameters (g=0.60, P=0.48) 
is shown in Fig. [4] (a), which includes the broadening and 
finite temperature effects. The calculated Fe 3d spectrum 
has the peak near Ep as a coherent quasi-particle peak, 
which becomes much sharper than in the PDOS with- 
out the electron correlation effect, and has two broad 
hump structures at around 2 and 4 eV. Thus, it repro- 
duces well the overall valence-band features experimen- 
tally obtained. The comparison between the experimen- 
tal and calculated Fe 3d states near Ep are shown in 
Fig. [5] The peak shift of the spectrum toward Ep due to 
the self-energy correction is obviously seen, the energy of 
which corresponds to KeS(uj) in the non-broadened spec- 
trum. One can also see that the leading edge behavior 
in the vicinity of Ep and the peak position is well repro- 
duced. We have measured the further higher-resolutioin 
spectrum at a low energy excitation (h^=40 eV) shown 
in Fig. [5] (b). In this photon energy the Fe 3d-state- 
dominant spectrum is obtained since the a of Fe 3d states 
is 16 times larger than that of Se Ap states^ Although 
the spectrum should be very sensitive to the surface elec- 
tronic structures, the spectral features agree well with 
the calculated Fe 3d spectrum. 

There is a discrepancy between the experimentally and 
theoretically obtained Fe 3d spectra in the Eb =0.2-3 eV 
as seen in Fig. |U[a) and FiglSJa). This should be due 
to the contribution of the incoherent part of the spec- 
tral function, A 1 d ac (k,Lu). In the present analysis, only 
the pole part of the Green's function, that is, the coher- 
ent part is considered. Thus, the quasi-particle spectral 
weight is reduced to Z kl 



+00 



dcuA d oh (k,ij) = Z k {< 1) 



(2) 



where Z k and A' d oh (k,u) are the renormahzation factor 
and the coherent part of the spectral function A d (k,oj) 
for Fe 3d states. According to the sum rules of the spec- 
tral function, the following spectral weight of the incoher- 
ent part remains in the occupied and unoccupied states, 



duAT 



l-Zu 



(3) 



The incoherent part should appear at hu> ~ e k +Rei7(w) 
on the higher Eb side of the quasi-particle peak, for in- 
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TABLE I: Experimentally and theoretically obtained mass enhancement factor (Z 1 or m* /rat) for typical Fe-based super- 
conductor systems, f indicates the non-superconductor. The Sommerfeld coefficient 7 is also listed. 



LaFePnO 

m*/m b [PES] 1.5-2.2 (Pn=P), f.8 (Pn=As, F-doped) 

Ref. [1^,14 ,23] 

m*/m b [DMFT] 1.9-2.2(Pn=P), 1.6 (Pn= f As) 

Ref. [40,41] 

7 (mJ/mol K 2 ) 10.1-12.5(Pn=P) 

Ref. [43,44] 



^BaFeaAsa (Ba,K)Fe 2 As 2 FeSe Fe(Se,Te) 

L5 2?7 3l5 6^20 

Ref.[13j Ref.[13j Ref.[3j 

1.8-2.1 2-5 

Ref.[42j Ref.[llJ 

16-37 23 5.4-9.1 39 

Ref.[45J6j Ref.[46] Ref.[3J] Ref.[16] 



stance, fiu ~ l.OeV for e , ~0.65 eV in the largest peak 
of the original Fe 3d PDOS. This is consistent with the 
energy region in which the disagreement between exper- 
imental and calculated spectra are seen. In fact, in very 
recent articles, theoretical works predict that the inco- 
herent spectral weight appears around Eb —2 eV as a 
lower Hubbard band . 10 ' 11 

The Z^ 1 at k = hp, that is, Z^ 1 depends on the real 
part of the self energy ReI7(u;) as follows,— 



2-1 sl _ dReS(u;) 



= 1 + 7*. 



(4) 



The Z~ x represents the mass enhancement due to the 
band narrowing of the quasi-particle peak. In the present 
work, we have obtained Z =3.6, which is about twice 
as heavy as the reported values in LaFeAsOo.94Fo.o6<^ 
In Table H the obtained Z -1 , that is, m* /nib (the ratio 
of the enhanced effective mass to the bare band mass) 
and some reported values related to the electron corre- 
lation in Fe-based superconductor systems are summa- 
rized. Spectroscopic and theoretical studies clearly in- 
dicates the common tendency of the electron correlation 
strength in these systems, that is, stronger on the right 
side, although the 7 value in FeSe is peculiarly smaller 
than in other superconductors. We conclude that the 
electron correlation in the FeSe superconductor is rela- 
tively strong, being quantitatively consistent with the re- 
sults of the DFT +DMFT (the density-functional theory 
combined with the dynamical mean-field theory) study.— 



IV. SUMMARY 

We have performed the bulk-sensitive SXPES and 
HAXPES for Fe-based superconductor FeSe. It is sug- 



gested in the Fe 2p core-level PES that the tetragonal 
FeSe has the Fe 3c? electrons with the highly coherent 
character unlike the non-superconducting FeSe with the 
hexagonal structure. In comparison to the results of band 
structure calculations, the Fe 3d band narrowing and its 
energy shift are revealed, which originates from the elec- 
tron correlation effect. Considering the self energy cor- 
rection, the larger mass enhancement (Z -1 ~3.6) than 
in other Fe-As superconductors are obtained. In addi- 
tion, the incoherent part of the quasi-particle spectrum 
are found and successfully separated. The obtained mass 
enhancement value and the energy position of the inco- 
herent part are consistent with the results of recent dy- 
namical mean- field calculations. 
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